Spin-waves in antiferromagnetic single crystal LiFePC>4 
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Spin-wave dispersions in the antiferromagnetic state of single crystal LiFePCU were determined by 
inelastic neutron scattering measurements. The dispersion curves measured from the (010) reflection 
along both a* and 6* reciprocal-space directions reflect the anisotropic coupling of the layered Fe 2+ 
(S = 2) spin-system. The spin-wave dispersion curves were theoretically modeled using linear 
spin-wave theory by including in the spin-Hamiltonian in-plane nearest- and next-nearest-neighbor 
interactions ( Ji and J2), inter-plane nearest-neighbor interactions (Jj_) and a single-ion anisotropy 
(D). A weak (010) magnetic peak was observed in elastic neutron scattering studies of the same 
crystal indicating that the ground state of the staggered iron moments is not along (010) direction, 
as previously reported from polycrystalline samples studies, but slightly rotated away from this axis. 

PACS numbers: 75.25,+z, 75.30.Ds, 75.50.Ee 



I. INTRODUCTION 

Lithium-orthophosphates LiMP0 4 (M = Mn, Fe, Co, 
Ni) have attracted a renewed interest in recent years, 
both for their relatively high lithium ionic-conductivity 
that can potentially be applied in rechargeable battery 
technology Q], and for their intriguing magnetic prop- 
erties, in particular, the strong magneto-electric (ME) 
effect they exhibit 0|- In this regard, of particular im- 
portance is LiFePGi, as it has already been tested as a 
high-potential cathode in secondary Li-ion rechargeable 
battery H 0, IS El Like other members of the lithium- 
orthophosphates, LiFeP04 is an insulator adopting the 
Puma space group In this structure, the Fe 2+ ion 

occupies the center of a slightly distorted MOq octahe- 
dron that shares oxygen anions with a PO4 tetrahedron 
forming a closely packed oxygen framework. The Fe 2+ 
ions (S = 2) form corrugated layers that are stacked 
along the [100] crystallographic axis, as shown in Fig. 
H^a). Nearest neighbors in the b-c plane are coupled 
magnetically by a relatively strong exchange interaction, 
J\ through an Fe-O-Fe oxygen-bond, whereas in-plane 
next-nearest-neighbors are coupled (J 2 ) via Fe-O-O-Fe 
@ (see Fig. □» for definit ion of exchange couplings). 
Interlayer magnetic coupling is mediated by a phosphate 
ion via an Fe-O-P-O-Fe bonding [To]. Thus, the olivine 
family of LiMPC>4 exhibits highly anisotropic proper- 
ties which are between those of two- (2D) and three- 
dimensional (3D) systems [Til [l2f . 

The magnetic properties of IAMPOa svstems have been 
studied since the early 1960s [ll Ej, Ell- They all un- 
dergo an antiferromagnetic transition at low tempera- 
tures to a similar magnetic arrangement differing only in 
the orientation of the staggered spins. Nearest- neigh- 
bor (NN) spins in the b-c plane are anti-parallel and 
the stacking along the a-axis is such that ferromagnetic 
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sheets perpendicular to the 6-axis are formed; nearest 
neighbor sheets are anti-parallel giving rise to (010) fun- 
damental magnetic reflection that, depending on the di- 
rection of the magnetic moment, can be intense or ex- 
tinct. Earlier neutron diffraction studies of polycrys- 
talline samples showed that the magnetic space group 
of LiFeP04 and LiCoPCU is Puma' with the spins ori- 
ented along b crystallographic direction (i.e., the (010) 
reflection is absent), and PnrrCa magnetic space group 
for LiNiP04 and LiMnP04 with the spins aligned paral- 
lel to the c-axis (i.e., strong intensity at the (010) reflec- 
tion) [HEIEIEIEI 

Il8l . Recent neutron diffraction 
measurements of single crystal LiCoP04 reported finite 
intensity at the (010) reflection, interpreted in terms of a 
ground state with a spin direction that is slightly rotated 
from the &-axis 12]. Weak ferromagnetism has been re- 
ported for LiMnPG-4 EH and LiNiP0 4 EH at a temper- 
ature below TV. 

Recently, the magnetic structure and properties of 
lithium orthophosphates have been reexamined theoret- 
ically and experimentally [2! |2l|. Rousse et al. [2l| 
reported neutron diffraction results from polycrystalline 
samples confirming the collinear structure below Tjv 
= 52 K. Magnetic properties of LiFeP04 investigated 
by Mossbauer spectroscopy and ma gne tization measure- 
ment determined that Tn — 50 K [22|. Theoretical es- 
timations of the spin exchange coupling by spin-dimer 
analysis, while neglecting the single ion-anisotropy, yield 
the following values J\ = 1.08 meV, J2 = -0.4 meV 
and J x = -0.92 meV @. However, so far there has 
been no experimental determination of the exchange cou- 
pling among Fe 2+ spins and of the single-ion anisotropy 
in LiFeP04, for comparison with the theoretical predic- 
tions. Knowledge of exchange couplings and single- ion- 
anisotropy is also important for understanding the origin 
of the strong magneto-electric (ME) effect in LiFeP04 
E3- All lithium orthophosphates exhibit a strong yet 
anomalous linear magnetoelectric (ME) effect with re- 
spect to the observed ME tensor components, a xy , a yx , 
for LiFeP04 and LiCoPC>4, and a xz , a zx for LiNiP04, 
as expected with their respective antiferromagnetic point 
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FIG. 1: (color online) (a) Atomic structure of LiFeP04. The 
Fe 2+ ions form buckled layers stacked perpendicular to the 
[100] crystallographic direction. The ground state of LiFePCU 
is collinear antiferromagnetic with the average moment along 
b direction, (b) Spin arrangement of the two Fe 2+ layers, 
the in-plane nearest and next-nearest neighbor interactions 
Ji and J2 and inter-plane nearest neighbor inaction Jj_ are 
labeled. 



groups mmm' and mm'm [23|, |2J, |2jJ 26] . In particular, 
the ME effect measurements of LiFeP04 as a function 
of temperature reveal a decrease of the ME coefficient 
alon g on e direction a yx (T) below a maximum close to 
Tn |23t Detailed determination of the magnetic 

structure using neutron diffraction from single crystals 
can shed light on the origin of these anomalies. 

Herein, we report measurements of spin-wave disper- 
sion curves of single crystal LiFePC>4 by inelastic neutron 
scattering measurements. Spin wave dispersion curves 
can provide the values of exchange interactions and the 
single ion anisotropy. The measured dispersion curves 



were modeled using linear spin wave theory by including 
the in-plane nearest-neighbors (NN) and next-nearest- 
neighbors (NNN) spin couplings the inter-plane nearest- 
neighbor spin coupling and the single-ion anisotropy. We 
have also employed single-crystal elastic neutron diffrac- 
tion techniques to investigate whether there are subtle 
deviations from the previously reported magnetic struc- 
ture determined from neutron diffraction measurements 
of polycrystalline samples. 



II. EXPERIMENTAL DETAILS 



LiFeP04 single crystals were grown by standard flux 
growth technique (LiCl was used as the flux) from a 
stoichiometric mixture of high purity FeCi2 (99.999% 
Aldrich) and Li 3 P0 4 (99.999% Aldrich) fzft. The grown 
single crystals have a dark greenish color. The compo- 
sition and structure were confirmed by carrying out Ri- 
etveld analysis of the X-ray powder diffraction (XRD) 
data, using the GSAS software package [2j|. No extra 
peaks from impurities were detected in the XRD pat- 
tern. Powder, for the XRD, was produced by crushing 
typical isolated single crystals from the melt. The lattice 
parameters yielded from the refinement (a = 10.337 A, 
b = 6.011 A, and c = 4.695 A) are in good agreement 
with literature values [i H I20 I l2l | . 

Neutron scattering measurements were carried out on 
the HB1A triple axis spectrometer at High Flux Isotope 
Reactor at Oak Ridge National Laboratory. A monochro- 
matic neutron beam of wavelength A = 2.366 A (14.6125 
meV, k a = 2vr/A = 2.656A" 1 ) was selected by a double 
monochromator system, using the (0 2) Bragg reflec- 
tion of highly oriented pyrolytic graphite (HOPG) crys- 
tals. HOPG crystal was also used as analyzer for both 
the elastic and the inelastic studies. The A/2 component 
in the beam was removed by two HOPG filters located 
before and after the second monochromator. The col- 
limating configuration 40'- 40' - Sample - 34 '- 68 ' was 
used throughout the experiments. Temperature measure- 
ments and control were achieved by a conducts LTC- 
20 using Lake Shore silicon-diode temperature sensors. 
An irregular shaped single crystal (weight ~ 0.4 g) was 
folded in aluminum foil and mounted on a thin aluminum 
post. It was then sealed in an aluminum can under he- 
lium atmosphere and loaded onto the tip of a closed-cycle 
helium refrigerator (Displex). Two temperature sensors 
were mounted on the cold-tip of the Displex and on the 
sample can. The temperature was controlled using the 
cold-tip sensor. The temperature difference between the 
two sensors was about 0.2 K over the temperature range 
investigated. 
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III. RESULTS AND DISCUSSION 
A. Elastic neutron scattering 

The LiFePC>4 crystal was oriented such that the a-b 
plane coincided with the horizontal scattering plane of 
the spectrometer. The elastic measurements confirmed 
that the magnetic structure of LiFePC>4 is antiferromag- 
netic with the main direction of the moment oriented 
along the 6-axis. However, contrary to the previous neu- 
tron diffraction experiments performed on powder sam- 
ples, we have detected the presence of the (0 1 0) reflec- 
tion. The intensity of this peak is relatively weak but 
its intensity follows a similar temperature dependence as 
that of a stronger magnetic peak (210) (see Fig. 0). ft 
is worth noting that the (0 1 0) peak is forbidden by the 
symmetry of the Prima! magnetic space group previously 
propos ed to describe the spin arrangement in LiFeP04 
|15l l2Cj . In the Pmna crystal symmetry, the point group 
of the Fe 4c site is m y and the only allowable magnetic 
point groups are, therefore, rriy with the Fe magnetic 
moments along the 6-axis (perpendicular to the mirror 
(010) plane) and m' v with the magnetic moment lying 
in the mirror (010) plane. The magnetic contribution to 
the (010) peak indicates that the ordered moment is not 
strictly oriented along the 6-axis and a small component 
perpendicular to this axis is present. This implies a low- 
ering of the symmetry of the magnetic space group where 
both magnetic components (along and perpendicular to 
6-axis) are allowed. From the intensity ratio of the two 
reflections 7(010)/ 7(210) at low temperatures, we can es- 
timate the angle of the staggered moment with respect 
to the 6-axis, by using the following relation, 



F(ow) sin(a(oio)) 



Ft 



(210) 



sin (a 



(210); 



I(ow) sin( 2 ^(2io)) /(2io) 
y ^(210) sin(20(oio)) /(oio) 
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where, i^io) an d F(oio) are the magnetic structure fac- 
tors of peaks (210) and (010), ck( 2 io), ck(oio) are the angles 
between the scattering vector of reflections (210), (010) 
and the ma gnet ic moment, and fi are the corresponding 
form factors |12j. Using /(2io)//(oio)~ 0.85, we estimate 
the moments are rotated by 7.5 ± 0.5 deg toward c-axis 
or 3 i 0.5 deg toward a-axis. The ratio between the 
magnetic and nuclear contributions to the peak intensi- 
ties of reflections (I ma g / Inuc) can be used to determine 
the average magnetic moment, /i from 
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where for the reflection in question, F nuc , and F mag are 
the nuclear and magnetic structure factors, I nuc and I ma g 
are the nuclear and magnetic intensities and f(Q) is 
the magnetic form factor of Fe 2+ at momentum trans- 
fer Q = 2kosm9. I mag can be calculated from the peak 
intensity difference at temperatures above and below Tjv. 
Using the peak intensities of (210) at 300 K and 10 K to 
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FIG. 2: Temperature dependence of the square root of the 
normalized integrated intensity at two reflection peaks (210) 
(a) and (010) (b). The transition temperature obtained from 
the fitting is Tjv = 50 ± 0.5K and the critical exponent (5 = 
0.27 ± 0.03. 



calculate the I mag and I nuc , and using /(210) = 0.85, the 
calculated average magnetic moment n for Fe is 3.93 ± 
0.05 /is, which is very close to the values of 3.99 hb and 
3.8 /is obtained in Refs. [2|| and |3fJ |. 

To determine the temperature dependence of the or- 
der parameter, the (010) and (210) reflections were mon- 
itored as a function of temperature in the range 10 to 
60 K. Figure Ufa) and (b) show the square root of the 
integrated intensity (%//) versus temperature. The \fl 
quantity is proportional to the antiferromagnetic stag- 
gered magnetization, i.e., the order parameter. It was 
fitted to the following power law function near the tran- 
sition temperature, 



(3) 



where, M is the sub-lattice magnetization at T = K, 
t = (1-T/Tn) is the reduced temperature, and (5 is the 
critical exponent. For the two magnetic peaks (010) and 
(210) the obtained transition temperatures are the same, 
= 50 ± 0.5 K and the critical exponent j3 is 0.27 ± 
0.03. The transition temperature is very close to the val- 
ues reported in the literature 50 ± 2 K [IE E3, E3 The 
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critical exponent f3 is slightly smaller than that calculated 
for the 3D Ising model (f3— 0.32) [31|. The temperature 
dependent background-like scattering above Tn and be- 
low ~ 60 K indicates some kind of critical scattering due 
to short-range order formation or due to a dimensionality 
cross-over. 



B. Inelastic neutron scattering 

The spin wave excitations were measured at 10 K along 
the (£, 1, 0) and (0, 0) reciprocal space direc- 

tions, for energy transfers (energy loss mode) ranging 
from 1 to 8.5 meV. As illustrated in Fig. E9(a,b,c), well 
defined dispersive magnetic-modes of resolution-limited 
energy-width were observed at all wave vectors. A typi- 
cal constant-q scan, performed at the zone center (010), 
is shown in Fig. y](a), indicating a single excitation at an 
energy transfer of 5.86 ± 0.04 meV. The increase in inten- 
sity at energy transfer below approximately 1 meV is due 
to the quasielastic scattering from the newly observed (0 

1 0) magnetic Bragg peak. Above the transition temper- 
ature (Tjv = 50 K), at approximately 55 K, the inelastic 
peak at 5.86 meV disappears, confirming its magnetic ori- 
gin. Such an energy gap in the dispersion curve is usually 
driven by single ion anisotropy. A similar energy gap of 

2 meV at 2 K, was also observed in the LiNiPC>4, and 
it was found to decrease with increasing the temperature 
|3^ . In the case of LiFePC>4, measurements performed 
at different temperatures indicate that the energy gap is 
temperature independent. The inelastic scattering signal 
measured at different constant wave-vectors £ along the 
(100) and (010) reciprocal-space directions, at 10 K, are 
shown in Figure (b, c). The data were fitted to Gaus- 
sian functions (solid line in Fig. [3] (a,b, c)) where the 
background was assumed to be constant. 

The spin- wave dispersion branches deduced from these 
fits, for both b* and a* reciprocal-space directions, are 
plotted in Figure^] It is shown that the dispersion curves 
montonically increase in energy with £, and that the spin- 
waves propagating in the plane along the (010) direction 
have higher frequencies than those propagating transver- 
sally along the (100). Qualitatively, this behavior reflects 
the anisotropy in the strength of exchange couplings in 
the system; as expected the in-plane exchange couplings 
are much stronger than those between planes. 

To construct a Hamiltonian for the spin system, we 
recall that in LiMP04 olivine family the in-plane super- 
exchange or super-super exchange interactions between 
nearest and next-nearest neighboring Fe 2+ ions (Ji and 
J2) are expected to be much stronger than that between 
the nearest inter-plane neighbors (J±) There- 
fore, we propose the following Hamiltonian, 
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FIG. 3: (color online) (a) Neutron scattering intensity as a 
function of the energy transfer E at (010) peak at 10 K and 
55 K. (b) Constant-Q scans taken at 10K, at different wave- 
vectors (0, 1+ £, 0) and (c) Constant-Q scans at (£, 1, 0). 



where, J% and J 2 are the in-plane NN and NNN cou- 
pling constants, respectively, and J± is the inter-plane 
NN coupling constant. The illustrations of Ji, J2 and J± 
are shown in Figure QJb). D is the single- ion anisotropy 
constant quantifying the tendency of the spins to align 
along the easy axis (the S z component is defined to be 
along the direction of the moment at the ground state 
- &*). The Ising like ground state of the system is be- 
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lieved to be invoked by the single-ion anisotropy term 
which comes about from crystal field effects and spin- 
orbit coupling [341. Using the antiferromagnetic spin- 
wave theorv [35t l36| . the lattice with JV sites was divided 
into two sublattices A and B, where nearest-neighbors of 
an Fe 2+ site in one sub-lattice are all sites in the other 
sublattice. The next-nearest neighbors of an Fe 2+ site 
are in the same sublattice. The magnon dispersion curves 
were calculated using the Holstein-Primakoff spin oper- 
ator transformation to linear approximation (i.e., linear 
spin- wave theory [13 )• The resulting spin- wave disper- 
sion is given by 



hw = A 2 - F 2 



(5) 



where A = (2J 1 ZS - 2J 2 ZS - 2J±ZS + 2J 2 ZS-/ 3N + 
2J ± ZSj± + 2DS) and F = 2J 1 ZS-f 2N , in which Z is the 
number of the nearest neighbors Z = 4, S =2 for Fe 2+ . 
72JV, 737V and j± are calculated using the following equa- 
tion: 



7(2jV,3jV,_L) 



,iQr 



(6) 



where r = (5, £, S±) are the components of vectors to the 
intra-plane nearest, next-nearest neighbors and to the 
inter-plane nearest-neighbor. We get 



72JV = cos(7rfcy) cos(7tA; 2 ) 



73JV = -^(cos(2nk y ) + cos(2nk z j) 



(7) 
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FIG. 4: (color online) Spin-wave dispersion curves along the 
6* and a* reciprocal space directions. Solid lines are fits ob- 
tained from Linear Spin- wave Theory using Eq(5). 



LiFePC>4. Our results show that although there are 
competing interactions between NN and NNN spins 
in LiFeP04, they do not lead to more complicated, 
incommensurate or non-colinear, magnetic structures. 
This is in contrast to the observation of incommensu- 
rate magnetic phases in LiNiP04|2^|. These compet- 
ing interactions may explain the observation of weak- 
ferromagnetism in LiMPC>4 svstems|l8j. They may also 
be related to the observation, in this study, that the stag- 
gered magnetic moment is not aligned along a princi- 
pal direction. From the gap in the spin wave disper- 
sion curve, we have been able to extract the single-ion 
anisotropy in LiFePC>4 using linear spin-wave theory. 



7_l = cos(wk x ) cos(irk y ) (9) 

The experimental data along the (0, 1+5, 0) and (6, 1, 
0) directions in Fig. 0] were fitted to Eq. (JSJ yielding 
the following values: J x = -0.662 ± 0.02 meV, J 2 = -0.27 
± 0.02 meV, J± = 0.021 ± 0.001 meV and D = -0.37 ± 
0.01 meV. 

The in-plane nearest- and next-nearest-neighbor cou- 
pling constants quantitatively agree with theoretical 
calculations 0, J\ = -1.08 meV and J 2 = -0.4 meV 0- 
The two spin couplings, J\ < and J 2 < 0, compete 
oppositely over the alignment of in-plane NNN spins; 
whereas J\ leads to parallel alignment of NNN spins J 2 
favors their antiparallel alignment. Such competing in- 
teractions can lead to incommensurate phases 38] which 
were not found in this system. However, incommen- 
surate phases have been reported for the isostructural 
LiNiP0 4 |2jj. The inter-plane coupling J ± = 0.021 meV 
determined in this study is significantly smaller than the 
theoretical one J ± = -0.92 meV 0. It should be noted 
that single-ion anisotropy was not considered in the the- 
oretical calculations |j- 

To summarize, we have measured spin- wave dis- 
persions and determined spin exchange couplings in 
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